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Abstract 

The ground state dynamics of an entropy barrier model proposed recently 
for describing relaxation of glassy systems, is considered. At late stages of 
evolution the dynamics can be described by a simple variant of the Ehrenfest 
urn model. Analytical expressions for the relaxation times from arbitrary 
initial states to the ground state are derived. Upper and lower bounds for the 
relaxation times as a function of system size are obtained. 
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Typeset using REVT^ 



The Ehrenfest urn model has played a crucial role in formulating and clarifying several 
fundamental and subtle concepts of statistical mechanics. In this model distinguishable balls 
are placed in two boxes (or urns). The dynamics consists of picking up a ball randomly and 
transferring it from its box to the other. Concepts like what does one mean by an equilibrium 
state, how does the system approach, eventually reach and thereafter persist in equilibrium, 
and the meaning of fluctuations that take the system away from its equilibrium state, become 
transparent when considering the time evolution in the Ehrenfest urn model. Interest in the 
Ehrenfest urn model, its variants and its generalizations has been revived recently following 
the work of Ritort [Q; see for some subsequent work on Ritort's and related models. 

Ritort's model essentially describes the relaxation of a nonequilibrium system. The 
ground state dynamics of Ritort's model at late stages of evolution is describable by a 
simple variant of the classical Ehrenfest urn model 0. In this letter we employ first passage 
time formulation and obtain analytical expressions for the relaxation times from an arbitrary 
initial state of (the variant of) the Ehrenfest urn model. 

In the model proposed by Ritort, distinguishable balls are placed in boxes. Energy 
is defined as minus the number of empty boxes. The dynamics is defined as follows. Select a 
ball and a box randomly and independently. If the selected box is non-empty deposit the ball 
in the box. If the box is empty, the transfer of the ball would increase the energy; hence carry 
out the transfer only with a probability given by the Boltzmann factor exp[— Aii^/i^'sT], 
where AE (= 1) is the increase in the energy, Kb is the Boltzmann constant and T is the 
temperature. The static properties of this model can be calculated exactly. For example at 
the ground state, which is A^- fold degenerate, the energy per particle is —1 + All the 
balls are in one of the A^ boxes at the ground state. Ritort's model could explain several 
important characteristics of the glassy systems, like anomalously slow relaxation, aging, 
and hysteresis. Ritort's model is perhaps the first model in which only entropy barriers 
are present explicitly. The ground state djTiamics of the model is equivalent to an A^ - 
urn generalization and variant of the Ehrenfest model. Moving a ball to an empty box 
is disallowed. Thus once a box becomes empty it stays empty for ever. As a result the 
number of empty boxes increases with time. The energy decreases with time and eventually 
reaches its ground state value. The relaxation of the energy is slow and becomes slower 
with time. This is easy to understand at least qualitatively, since lowering the occupancy 
of a box becomes less and less likely as the number of balls in the box becomes fewer and 
fewer. Let us denote by fl{N, N) the set of all possible configurations of distributing A^ 
balls in A^ boxes. Let r denote the time it takes for the system to relax to the ground state, 
averaged over the ensemble Q{N, N) of initial configurations. It is clear that the system 
before reaching the ground state, arrives invariably at one of the simpler non equilibrium 
configurations with only two boxes non empty. The other N — 2 boxes have since been 
emptied during the evolution. Let us denote by Q{2, N) the set of all possible states with N 
balls in 2 boxes. Also let ri denote the relaxation time averaged over the ensemble ^2(2, N) 
of initial states. In a recent work, Lipowski |^ showed that Ti is less than but almost equal 
to r, for large A^. 

The purpose of this letter can be stated as follows. Lipowski p derived analytical 
expression for the relaxation from one of the simple states belonging to the set Q{2,N). 
This state consists of only one ball in one box and the remaining A^ — 1 balls in the other. 
Relaxation times from other states like k balls in one box and the rest in the other, were 



2 



calculated through a set recursion relations. Lipowski's work is a clever application of the 
formulation of Darling and Siegert |^ , for first passage times of random walks. In this letter 
we employ a simpler and more recent first passage time formulation |^ . This method enables 
us to obtain closed-form expression for the relaxation time from an arbitrary initial state, a 
task which has not been possible in Lipowski's formulation. More importantly, employing 
the analytical expressions, we establish useful upper and lower bounds for the relaxation 
times from various initial states in terms of the system size. 

We first consider the case with even number of balls distributed in two boxes. Accordingly 
let 2N be the total number of balls in the system. Let us consider a configuration with k 
balls in one box and the remaining 2N — k balls in the other. We say the system is in state 
k if the number of balls in the box of lower occupancy equals k. Thus k is less than or equal 
to A^. Let Gk,k-i{^) denote the probability for the system to go from the state k to the state 
k — 1, for the first time in exactly u time steps. Thus u is a discrete random variable and is 
called the First Passage Time(FPT). The first passage from k to k — 1 can happen in two 
ways: (a) Select a ball from the box containing k balls and transfer it to the other taking 
one time step. The probability for this is k/2N. At the end of this step we are in the target 
state k — 1. Hence u, the FPT is unity, (b) Select a ball from the box containing 2N — k 
balls and move it to the other box taking one step. The probability for this is {2N — k)/2N. 
At the end of this step we are in the state k + 1. Now make a first passage from the state 
+ 1 to the state /c — 1 in the remaining z/ — 1 steps. The above considerations hold good for 
all A; < A^ — 1. For k = N, we find that since both the boxes contain A^ balls each, selecting 
a ball from either and moving it to the other leads to the target state k — 1, and the FPT 
is unity. The state A^ is thus reflecting. We have therefore, 

k 2N — k " 

Gk,k~i{v) = i^^v,! + Gk+i,k-i{,v - 1), for k= 1,2, ■■■,N-1, 

Gn,n^i{i') = (1) 

The above is the complete set of A^ equations for the FPT densities. Since the system starting 
from state k + 1 can not reach the state k — 1 without visiting the state k, Gk+i,k-i{i^) in 
the above can be expressed as a convolution given by, 

Gk+i,k-iii^) = Gk+i,k{il)Gk,k-i{v - v)- (2) 

r?=l 

The convolution in the above holds good in general for one dimensional problems with 
nearest neighbour hopping. 

When the number of balls is odd, say 2N — 1, the equations for the FPT densities 
remain formally the same except that, now the state A^ — 1 is reflecting, in the following 
sense. Consider the first passage from the state A^ — 1 to the state N — 2, in u steps. There 
are two ways: (a) Move a ball from the box containing A^ — 1 balls to the other box, in one 
step; the probability for this is (A^ — 1)/(2A^ — 1). At the end of this step we are in the 
target state N — 2. Hence the FPT is unity, (b) Move a ball from the box containing A^ 
balls to the other box, taking one time step. The probability for this is N/ {2N — 1). At the 
end of this step we are in state A^ — 1, the same state we started with. Now make a first 
passage from this state A^ — 1 to the target state A^ — 2 in the remaining u — 1 steps. Thus 
for odd number of balls we get. 
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k 2N — 1 — - 

gfc,fc-i(z/) = ^^^^.,i+ Gfc+i.fc-ilt^-l), for k=l,2,---,N-2 

iV - 1 

G'^_i,7v-2('^) = 2jv_ i ^'^.i + 2iV- - (3) 

The above constitute the complete set of iV — 1 equations for the FPT densities for the case 
with odd number of balls. To solve the recursion relations and (|^) for the first passage 
time densities we employ generating function technique. 

Let Gij{Z) denote the generating function for the FPT, defined as, 

oo 

G,,(Z) = 5:Z^G,,(z/). (4) 

u=l 

Multiplying both sides of Eqns.(|l]) by Z'^ and summing over u from 1 to oo, we get, 

k 2N — k 

Gk,k-i{Z) = Z— + Z^^Gk+i,k-i{Z), for k=l,---,N-l 
Gn,n-i{Z) = Z. (5) 

Equivalent relations, not given here, can be obtained for the case with odd number 
balls. From Eqns.(^, a terminating continued fraction relation for Gk,k-i{Z) can be derived 
by noting that Gk+i,k-i{Z) = Gk+i,k{Z) x Gk,k-i{Z), by virtue of convolution theorem. 
Substituting the convolution in Eqns. @, we get, 

Gk,k-i{Z) = J_P — , for k=l,..-,N-l 

Gn,n-i{Z) = Z. (6) 

In fact, by convolution we have Gmfi{Z) = JlfcLi Gk,k-i{Z), for m = 1, ■ ■ ■ , N. Thus in prin- 
ciple we have obtained the distribution of relaxation time from an arbitrary state belonging 
to Q{2, N), to the zero ground state, though the expressions are in Z space. 

For calculating the mean first passage time(MFPT), from the state k to the state k — 1, 
we differentiate Gk,k-i{Z) with respect to Z and set Z = 1. Let Fk^k-i denote the MFPT 
from k to k — 1. For the problem with even number of balls, we get, 

Fk,k-i = ^^fc ^ Fk+i.k + for k = 1, ■ ■ ■ , N - 1 

The above can be cast in a convenient matrix notation, 

\F) = A\F) + \U), (8) 

where \F) is a column vector (Fi_o -^2,1 ■ ■ ■ ,-^Af,Af-i)^ and \U) is the column vector rep- 
resenting the inhomogeneities, {2N/1 2N/2 ■ ■ ■ 1)^ Here the superscript ^ denotes the 
transpose operation. The N x N matrix A has elements given by Aij = x (2A^ — 
We can cast Eq. (|]) as B\F) = \U) where B = I — A. The matrix B has all its diagonal 
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elements unity; all the other elements except those in the first upper diagonal are zero. The 
matrix elements of B are given by, 



X 6,^-1 + Si 4. 



(9) 



We give below the matrix B explicitly, to facilitate easy visualization of solutions we are 
going to derive shortly. 



B 



( 1 






V 



2N-1 

1 





/ 2N-2 
V 2 



























-, _ ( 2N-{N-1) \ 
^ \ 2N ) 







(10) 



To calculate the m— th element Fm,m-i of the vector \F), we replace the m-th column of the 
matrix B by the vector \U). Let the matrix thus formed be denoted by S(m)- Then Fm,m-i 
is given by Cramer's rule. 



m,m— 1 



Tn)\ 



D[B] 



(11) 



where D[-] denotes the determinant. First we observe that the determinant of the matrix B 
is unity. The problem reduces to calculating the determinant of B(^rn) ■ 

Let us consider the case with m — 1. The determinant of can be easily written down 
by inspection as. 



-^1,0 = D 

+ 



2N 2N/2N-1\ 2N /2N -1\ /2N -2 



1 + 2 



+ 



1 2N /2N - 1\ f2N -2 
^2 ~N [ ~ 



'2N-{N-1) 
N-1 



(12) 



which can be cast in a compact form as sum over products given by. 



n=l 



n 



1 — -^^n,N 



n-1 

n 

k=l 



2N-k 



k 



N-1 ... 

n=l 



(13) 



where ^■''C(.) are the usual binomial coefficients. Noting that '^^Cn is the same as '^^C2N-m 
we see immediately, that 



2JV-1 



(14) 



n=l 
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Add the binomial coefficients ^^Cq = 1 and ^^C2n = 1 to both sides of the equation above. 
We find that the right hand side becomes 2^^, and we get, 



F,,o = 2^^-^ - 1, 



(15) 



which is precisely the expression derived by Lipowski ^j. Note that in Lipowski's paper N 
denotes the total number of balls whereas here the total number of balls is taken as even 
(2A^), see eq. (|l^) or odd (2A^— 1), see Eq. (|^). We have considered the two cases separately 
for bringing out clearly the subtle difference in the refiecting boundary while deriving the 
master equations, see Eq. (|l|) and Eq. (|^). 

Let us now derive closed- form expressions for -Fm,™,-!- To this end, we replace the m-th 
column of the matrix B by the vector | U) and construct the matrix -B(m) , whose determinant 
gives, 



N-l 



/2N 



m,m—l 



n-l 



2N-k 



k=m 



k 



1\ /2N 



N-l 



2N-k 



k=m 



k 



(16) 



First we multiply both sides of the above equation by 11^=1 (2 -k)/k= {m/2N) ^^C^, 
and get. 



m ^ m,m—l 



N-l 

E 

n=m 
N-l 



/'2N\''-^ 2N - k n \ /2N\^-^2N -k 



n 



n 

k=l 
1 

2 



2/ V iV 



n 

k=l 



2N r< 
Oat. 



(17) 



Now add to both sides of the above equation the term Y^n=i ^^C*„, and get, 

N-l 



\2A^ / 



n=l 



m-1 — E "^^^n + [ 

n=l ^'^ 



2N 



We see immediately that 



2 X 



'm— 1 



n=l 



2N-1 

E 

n=l 



2N 



(18) 



(19) 



If we add now ^^Cq = 1 and '^^C2n = 1, to both sides of the above equation, we find the 
right hand side is simply 2^^. We get 



m,m—l 



I2N 



m 



2Nn 

On 



m—1 



(2^^-^ - 1) - 5: 



n=l 



(20) 



It is easily seen that if we substitute m = 1 in the above we recover Lipowski's result P|, 
also derived explicitly in this paper, see Eq.(|T5| ). The relaxation time from any state k to 
the zero ground state can be obtained by summing Eq. (pOD over m from 1 to k. Thus we 
get. 



fc,0 



E 



f2N 



m—1 



(^22N-1 - l) _ ^ ^^Cr. 



n=l 



(21) 
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For odd number of balls the derivation proceeds in the same way, and we get 



fc,0 



E 

m=l 



2N -I 



m 



2N-1 



Cn 



(2 



2N-2 



1 



m— 1 

E 

n=l 



2N~l 



(22) 



Now that we have analytical expression for the relaxation time, Fkfl from an arbitrary 
state, we can estimate how much it deviates from Fi q, when the system size goes to infinity. 
From Eq. (EOl) it is clear that 



-^1,0 > -^2,1 > -^3,2 > ■ ■ ■ > Fn^i^n~2 > Fn,N-1 

In fact for large A^, we have, from Eq. (|20|), 



1). 



m,m—l JV- 



m—l,m—2 



which implies that. 



Fm,m-1 N^oo Fifl 



(m 



m — 1 
2N 



IV. 



(2iV) 



m—l ' 



(23) 



(24) 



(25) 



Since, Fi 



kfl 



= Fifl + ^2,1 H h we have, 

Fkfl mZoo Fifl 



— ^ 

^ {2Ny 



+ 



(2A^)^ 



+ ■■■ + 



(2iV)'=-i 



(26) 



for all k. Thus to the order of ^, we have 

Fkfl AtZoo -^1,0 



1 + 



2NJ 



for all > 2 and the correction is independent of fc. In fact it is easily shown from Eq. 
that, 



(27) 

m 



Fifl 1 



1 

2N 



< Fkfl <Fi,o 1 



1 

N 



(28) 



for all k > 2, when is large. Thus it is clear that, indeed Fifl is the principal time scale 
in the problem and relaxation time from any other state is only negligibly greater than Fifl 
for large systems. 

This letter, in a way complements the work of Lipowski f^. We have obtained analytical 
expressions for the relaxation times to the ground state (with all the balls in one box and 
the other box empty), starting from an arbitrary initial state (with say, k balls in one box 
and the rest in the other). We have shown that the relaxation time from the simple state 
k = 1 (with one ball in one box and the rest in the other) sets the principal time scale in 
the problem; relaxation from other states takes negligibly more time than this, for large 
systems. A natural question that arises in this context is whether we can construct a simple 
two-urn analogue of Ritort's model. To this end we need to suitably modify the energy 
function defined over the states of f2(2, A^). For example we can define energy as minus the 
absoluted value of the difference of the number of balls in the two boxes, i.e., Ein] = 2n — N, 
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where n < N/2 is the number of balls in the lower occupancy box(defining the state of the 
system) and is the total number of balls. It is easily seen that the state n = is the 
minimum energy ground state with £^(0) = —N. The maximum energy state is n = N/2, 
with E{N/2) = 0. It can be shown that an arbitray state k relaxes to the ground state in 
a time given by NYlm=i a-t zero temperature. However as the temperature increases, 
the relaxation of energy to its equilibrium value (at that temperature) becomes faster. It is 
indeed worthwhile investigating if other features like hysteresis and aging are also present 
in this simple model. 

KPNM thanks Forschungszentrum Jiilich, for the hospitality extended to him at the 
Institut fiir Festkorperforschung. 
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